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Abstract. The phi meson production on hyperons, ttY — * <f>Y and anti-kaons 
KN — ► <fiY is argued to be a new efficient source of phi mesons in a nucleus- nucleus 
' collision. These reactions are not suppressed according to Okubo-Zweig-Izuka rule 

in contrast to the processes with non-strange particles in the entrance channels, ttB 
and BB with B = N, A. A rough estimate of the cross sections within a simple 
hadronic model shows that the cross sections of ttY — > cj)Y and KN — > <fiY reactions 
can exceed that of the ttN — * <fiN reaction by factors 50 and 60, respectively. In the 
^ ■ hadrochemical model for nucleus- nucleus collisions at SIS and lower AGS energies we 

£h . calculate the evolution of strange particle populations and phi meson production rate 

due to the new processes. It is found that the catalytic reactions can be operative if 
the maximal temperature in nucleus- nucleus collisions is larger than 130 MeV and the 
collision time is larger than 10 fm. A possible influence of the catalytic reactions on the 
centrality dependence of the <f> yield at AGS energies and the (/) rapidity distributions 



(N 



(T) . at SPS energies is discussed. 

o 

ON 

o 



PACS numbers: 25.75.-q,25.75.Dw,25.80.Nv,25.80.Pw 
1. Introduction 



The (j) meson production is an important part of the study program of different nucleus- 
nucleus collision experiments in the whole range of collision energies: at the AGS [1,2], 
the SPS [3-7], RHIC [8,9], and even at deeply subthreshold energies at the SIS [10]. As 
the longest- living vector meson, the <fi is considered to be a good probe of the collision 
dynamics: it would decay mainly outside the fireball and the daughter hadrons would 
be weakly affected by rescattering. The dominant hadronic decay is <p —>■ K K, hence, 
the production of 0's interplays with the production of strange mesons in nucleus- 
nucleus collisions and can provide a complementary information on the strangeness 
dynamics. As other neutral vector mesons, the <ft contributes to dilepton production 
via decays <p — > e + e~ and fi + fi~ . The comparison of <fi meson yields in the hadronic 
and electromagnetic channels [11-13] can reveal interesting information about the later 
stage of the collision [14-17]. 

According to the SU(3) quark model the dominant component constituting the 
meson is a spin-one bound state of s and s quarks. Hence, the hadronic interactions 
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of phi mesons are subject to Okubo-Zweig-Iizuka rule, which tells that the interactions 
of a pure (ss) state with non-strange hadrons are suppressed. In practice, the strict 
implementation of this rule would imply vanishing of NN and (fiirp couplings. Indeed, 
the OZI-forbidden reactions are typically orders of magnitude smaller than the OZI- 
allowed, e.g., the ratio of uj to meson production cross section is about [18] 

a{7iN^<j)N) y ' 

This value is, however, three times smaller than the one which follows from the SU(3) 
symmetry and the experimentally known deviation from the ideal mixing between singlet 
and octet vector mesons [19]. The phenomenological models are able to reproduce the 
experimental cross section of the ttN — > <pN reaction by taking into account the OZI- 
violating (fipn coupling adjusted to reproduce the — > pir decay width [20,21]. The role 
of baryon resonances in this process was studied in ref. [22,23]. 

Since the OZI suppression weakens the production only by the ordinary hadronic 
matter and is lifted in the quark-gluon medium, the strong, of the order indicated in flTJ), 
enhancement of the yield was proposed in [24] as a signal of the quark-gluon plasma 
formation. An enhancement of the yield was indeed observed in experiments albeit 
to a lesser degree: the ratio of the yield in a nucleus-nucleus collision normalized by 
the number of participating nucleon pairs to the yield in a proton-proton collision lies 
between 3 and 4 for AGS/SPS energies [2,7]. Within the hadrochemical model [25] the 
factor 3 of enhancement in production could be explained only by assuming a decrease 
of hadron masses in hot and dense hadronic matter. In [25] the main contribution 
to the yields are given by the OZI-allowed process with strangeness coalescence 
KK — ► 0p and KA — ► 0iV . The string-hadronic transport model UrQMD [26] , in 
which the is produced mainly via KK — > reactions and no in-medium effects are 
included, systematically underestimates the production at SPS energies, cf. figure 9 
in ref. [7]. The experiments suggest some similarity in production mechanism of 
mesons and kaons. In Au+Au collisions at 11.7 A GeV/c [2] and in Pb+Pb collisions 
at 158A GeV/c [3,4] it was observed that the yield grows with the centrality more 
strongly than linearly in a similar way as the K + and K~ yields do. In contrast, the 
p + u) production scales linearly with a number of participating nucleons [4]. 

At the SIS the production was studied at the beam energies about 2 GeV per 
nucleon, which are less than the threshold energy of production in the nucleon-nucleon 
collision 2.6 GeV. The FOPI collaboration measured 0s in Ni+Ni reaction at 1.93 
AGeV [10]. The reported ratio of to K~ multiplicites is (<f>)/(K~) « (0.44±0.16±0.22) 
or (1.7 ± 0.6 ± 0.85) depending on the temperature of the emitting source (130 MeV 
and 70 MeV, respectively) which has to be assumed for an extrapolation to a full 
solid angle. These results are considerably larger (at least by factor 4) than the very 
preliminary estimates announced in [27]. Such a large abundance cannot be described 
by the transport model [28] where 0s are produced in reactions BB — > BBcj) and 
7rB — > (pB (B = N, A) with the dominant contribution from pion-nucleon reactions. 
The calculations with the IQMD code [29] underestimate the production data [10], 
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Table 1. Relative multiplicities of strange particles observed in nucleus collisions 
at different colliding energies and the estimate of the lowest values of irY — > <\>Y 
and KN — ^ <f>N cross sections at which this <p meson production mechanism could 
dominate the irN — > (f>N reaction (a(wN — > 4>N) ~ 0.03 mb). The empirical date are 
taken from [33] . 



Kab 




2 GeV/A 


4 GeV/,4 


6 GeV/A 


8 GeV/A 


10 GeV/A 


(K+) 

Wpart 

(A) 

Wpart 

2{K-) 

(TT++7T-) 




1.1 x 10~ 3 
1.3 x 10~ 3 


7.0 x 10~ 3 
8.4 x 10~ 3 
5.8 x 10~ 3 


14 x 10~ 3 

15 x 10~ 3 
14 x 10" 3 


23 x 10~ 3 
20 x 10~ 3 
22 x 10~ 3 


32 x 10~ 3 
30 x 10~ 3 
34 x 10~ 3 


<t(ttA — 
cr(KN - 


- 4>Y) > 
-> cf>Y) > 


23 mb 


3.5 mb 
5.2 mb 


2 mb 
2 mb 


1.5 mb 
1.3 mb 


1 mb 
0.9 mb 



too. Even the inclusion of other resonance channels pB —>■ <pN and 7riV(1520) — » (f>N 
studied in ref. [30] would not fully accommodate the results [10] on (p multiplicites. Note 
that the strangeness coalescence process could not contribute much to the yield at 
these energies since kaons have a long mean free path and most likely leave the fireball 
right after its creation without any further interaction. 

The large value of (4>)/{K~) reported in [10] seems to be confirmed by further 
analyses of the data [31]. The preliminary results by HADES experiments with Ar+KCl 
collision at 1.756 AGeV give also a comparable value for {$) / (K~) ~ 0.37 ± 0.13 [32]. 

In this paper we would like to address some puzzles of the <p production in heavy-ion 
collisions at various energies and propose a new type of the (f) production mechanism — 
the catalytic (ft production by strange particles, e.g. the reactions 

nY -> <j)Y , KN ^(j)Y, Y = A , E . (2) 

In contrast to the strangeness coalescence reaction, here the strangeness does not hide 
inside the <fi but stays in the system and the presence of K mesons is unnecessary. The 
efficiency of these reactions should be compared with the process tcN —>■ <f)N, which is 
found to be dominating in [28]. The reactions (T5]) are OZI allowed, so we win in cross 
sections compared to ttN — > (f>N a factor of the order of 75, cf. ([1]). We lose, however, in 
the smaller concentration of hyperons and anti-kaons compared to nucleons and pions, 
respectively. 

There is no experimental information about the cross sections for the reactions (T5]), 
therefore we will have to rely on some modeling. Let us first estimate how large the cross 
sections for these reactions should be so that their contributions become comparable to 
ttN — > (f)N7 This can be deduced from the following relations 

a(7rA -> (j)Y) > 6t(vtA -> <pN) , 

a(KN - <j>Y) > fr^±l ) a ^N - <j>N) . 
Here iV part number of nucleons participating in the collision, (A) , (K~) , (ir + ) , and 
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7T,q <f>, q tt, q </>,<? 




Figure 1. Diagrams contributing to a ttY — > 4>Y reaction 

(7r~) are the multiplicities of the corresponding particles produced in the collisions. 
Using mid-rapidity data for collisions at AGS energies, cf. [33], listed in table [Hand the 
experimental value o~(itN — > <pN) ~ 0.03 mb we get an estimate for the lower bound 
of the 7rA — > cf)Y and KN — > cf)Y cross sections. We see that a few mb is enough to 
make the cathalitic meson production comparable with the conventional mechanisms 
at AGS energies. For SPS energies these bounds could be even lower, since there are 
more kaons and hyperons. 

Although encouraging these rough estimates are not convincing yet. Indeed, the 
important is not only how much strangeness is produced in the collision but also when 
it is produced and at what temperature. The catalytic reactions ([2]) will be operative 
if there is enough strangeness produced and the medium is still hot enough. To study 
various scenarios quantitatively, in section 2 we estimate the cross section of the reactions 
([2]) within a simple model. In section 3 we set up a hadrochemical model for strangeness 
production. The production is discussed in section 4. 

2. Estimate of nY — > <\>Y and KN — > <pY cross sections 

We consider the Lagrangian of lightest strange particles, i.e. hyperons and kaons, 
interacting with pions and nucleons and include the meson as a gauge boson of the 
local U(l) symmetry associated with the strangeness conservation 

C = A(i 9-7 + 0- 7- m A )A + £ a (z<9-7 + 0- 7- ™ E )£ a 
+ ±(9-2 0)^ {d + t<f)fK- 

+ CVas (A 7m 75 S a + £ a 7m 7 5 A) d^ a + CVse % e abc (^ a l, l^ b )dV 

+ CW (N 7m 75 r a S a (0" + i g ^)K+ {&* - i g$ 0") ^ S a r a 7 „ 7s N ) 

+ C KNA (N 7m 75 A(d»+i g,,, r) K+{d»-i g<t> r) A llx 7 5 N ) , (3) 

where S a and 7r a denote the isospin triplets of £ baryons and pions; kaon and nucleon 
isospin doublets are K = (K + ,K°) and N = (p,n). The coupling <pKK is fixed from 
the — > A A" decay, ^ = 4.64. The couplings CVae = 0.63 m" 1 , C^ee = — 0.56m 7 ^ 1 , 
Crna = -OJm^ 1 and Cknt. = 0.4m" 1 we take from the analysis of meson baryon 
interactions [34]. 
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2.1. Ti Y 



Y reaction 



In figure [T] the lowest-order diagrams contributing to the n Y - 
The corresponding amplitude reads 

T = 9<p C nYY i u Y (p; s)YS Y (P+ q) (q • 7) IhUy (p; s) e* (g; A) 
+ 9^C nYY iu Y (p; s) (g-7) 75<Sy(p-?) l^u Y {p\ s) e* (g; A) 
= 0* CWy %(p; s) J YY (q) u Y (p; s) e* (q; A) , 



K process are shown. 



S Y (p) 



p-^ + M 



s y (p) 



p-j + M 



(4) 



p 2 -M 2 ' ~ rvv p 2 -M 2 
Here we use the following notations: p and g are the 4-momenta of a hyperon and a 
pion, respectively, in the initial state, p and q are the 4-momenta of a hyperon and a phi 
meson in the final state; M and M are masses of the initial and final hyperons; u Y (p, s) 
denotes the Dirac bispinor of a hyperon Y with the momentum p and the spin projection 
s; e* (g; A) is the wave function of an outgoing meson with the polarization A. For the 
moment we do not specify in (j3J) the isospin coefficient for the particular reaction. Taking 
into account the on-shell relations for hyperons (p-j — M) u(p, s) = = u(p, s) (p-j—M) 
we can reduce the current J YY (q) introduced in (j4]) to the on-shell equivalent form 

.(g) = i (M + M) S Y (p + g) 75 + 7s S Y (p - q) Y} ■ 



j yy \ 



J£ Y (q) 7o = J YY (-q) and Jfi(q) g M = . 



The current has the convenient properties 70 
The cross section is determined by the amplitude squared, summed over the spin and 
polarization of the final hyperon (s) and phi meson (A), and averaged over the spin of 
the initial hyperon (s) 



s,s,X 

- a 2 r 2 - 
n 2 r 2 - 

9d> U ttYY 



\ Tr {(p • 7 + M) J^ YY {q) (p • 7 + M) J YY (-q)} 
(M + M) 2 [f YY (p, q,p) + g YY (p, q,p) + g YY (p, -q,p)] 



where 

fyyip, q, P) = ~ Tr {(p • 7 + M) ^ S Y (p + g) 75 (p • 7 + M) ^ S Y (p - q) 75} 
+ lTr{(p-7 + M)7 5 5y(p-g)7 M (p-7 + M)7 5> Sy(p + g)7 M } 

= {{ p + qy-M 2 ){{p-q) 2 -M 2 ) [MM(2p.p + ml) 
+ (M - M)(Mp-q + Mp ■ q) - 2p • (p - q)p • (p + g)l , 

9yy(P, q, P) = ~ Tr {(p • 7 + M) 7 M Sy(p + ?) 75 (p • 7 + M) 75 S Y {p + g) 7^} 

= (( P+ ^-M^ [ 2( ^^ A?2)(p - ?+A?M) 
- (p-p-2MM) (2M 2 + ml 
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°yy( s ) = t^t" — \ [ dx jF-p y (s, u(s, x)) (5) 



Using the relations 2 p • q = s — M 2 — m 2 , 2 p ■ q = M 2 + m\ — u, 2 p ■ p = M 2 + M 2 — t , 
and u + s + 1 = M 2 + M 2 + m 2 + m| we can express T YY as a function of two invariants 
s and u . Then the total cross section can be evaluated as (cf. [35]) 

32 7T S p n Y(s 

— _L 

where 

u(s, x) = \{ m l + ml + M 2 + M 2 -s) + ( M2 ~ m lH M2 ~ m l) _ 2pMs) ^_ (s) x 

pIy(s) = M s M2 f M2 ) « M s ) = i(( s ~ m l~ M2 ') 2 ~ Am l m 2 ) ■ 

In calculations of the particle production in heavy-ion collisions one usually uses 
isosp in- averaged cross sections which can be defined as the sum of cross sections for 
active (allowed) isospin channels divided by the total number of isospin combinations 
of colliding particles. In ttA — > channel the following three reactions are possible 
7r ±,0 A — > (fiT.^' . Since all these processes have the same cross section with the isospin 
coefficient equal to one, the isospin-averaged cross section is 

<7 v A-*4>e{s) = cr AS (s) . (6) 

In the 7rS — > channel we have six possible reactions 7r ± S =F — > </>X° , tc^H — > </>S ± , 
and — > . The entrance channel, 7r£, counts 9 possible isospin combinations. 
The isospin averaged cross section is then 

2 2 

^2-^2 (s) = g 0-tt+E— >4>X°(s) = ~ (Tss(s) . (7) 

The 7rS — > 0A channels is represented by three reactions tt^II 11 — > 0A and 7r°S° — > </>A 
out of nine available for the 7rE entrance channel. 

0V2-^a( s ) = - ^+£--^0) = ^sa(s) • (8) 

The isospin-averaged cross sections (l6fTf8"]) are shown in figure [2] as functions of the 
collision energy above the threshold As 1 / 2 = s 1 ^ 2 — s\(^ . For comparison, we depict the 
cross section for the n~p — > (f>n reaction in the parameterisation [29] (solid line) and [36] 
(dashed lines) scaled up by factor 5. We see that the cross sections of the phi production 
on hyperons are in average about 50 times larger than the production on nucleons (at 
As 1 / 2 ~ 100 Mevjj]. 

2.2. KN — > <pY reaction 

The diagrams contributing to KN — > (f>Y processes in the lowest order are depicted in 
figure [31 The reaction amplitude is 

| Note that the isospin average cross section of nN — > <f>N reaction is a(nN — > (f>N) = \ <r(n~p — > (j>n). 




Figure 2. Isospin averaged cross sections of the (j> meson production in irY — > <\>Y (left 
panel) and in KN — » (right panel) reactions. On the left panel two lowest lines 
depict the parameterisations of the 7r~p — > 0n reaction from ref. [29] (solid line) and 
ref. [36] (dashed line) in comparison with the experimental data from [37,38]. Both 
the paremeterizations and the data are scaled up by factor 5. 



K,q </>,q K,q <j>, q K, q (/), q 




Figure 3. Diagrams contributing to a KN — > 4> Y reaction 



i Tyn = 9<fi C KNY u Y (p] s) (7" - (q - q) ■ 7 D K (q -q)(2q- g) M ) 75 u N (p; s) e* (g; A) 
- 9<t> C kny u Y (P; s) Y S Y (P + q)(q- 7) 7s u N (p; s) e* (g; A) 
= 9<t> C KNY u Y (p; s) J YN (q, q) u N (p; s) e* (g; A) , 

= (<? 2 - m^)- 1 . (9) 

As before, we do not write explicitly the isospin factors for each reaction. They will be 
restored later. Using the on-shell conditions we rewrite the current J YN in the on-shell 
equivalent form 

4 N (q, q) — (M + M) (D K {q -q)(2q- q)» + 7" S Y (p + q)) 75 • 
The squared and spin summed amplitude is then 

Fyn = \ E T YN T YN = ^Tr {(p ■ 1 + M) l5 (p ■ + M N ) l5 } D 2 K (q - q) (2 q - qf 

s,s,\ 

+ lTr{(p-7 + M)7 5 (p-7 + M^)7 5 Sy(p-g) 7 4 D K (q - q) (2q-q)» 
+ -Tr{(p-'y + M)>y ll SY(p-q)7 6 (p-7 + M N )y B } D K {q-q) (2g-gf 
+ -Tr{(p-7 + M)7»Sy(p + q)75 (p • 7 + M N ) 75 S Y (p + q) 7^} 
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-a- (vav) 2 (2<1-V) 2 (P-P-MM N ) 

Hq-Q) *-™>k) 

({p-p-M M N ) (2 q - q) ■ (2p + q) + 2 (p • q) (p ■ q) - 2 (p • q) (p ■ q) 
~ 4 ((q - q) 2 - m\) ((p + q) 2 - M 2 ) 

The function g-p^ follows from g y y after the replacement M — > Mn • Using the 
kinematic relations 2p- q = s — M 2 — m 2 , 2 p ■ q = s — — m 2 K , 2 p • p = M 2 + — t , 
2 q ■ q = m 2 ^ — m 2 — t,2p-q = s + t — Mfj — m 2 ^ , and 2 p • q = s + t — M 2 — m 2 we 
can express Tyn as a function of invariants s and t . The cross section of KN — > Y(f> is 
given by 

<Tyn(s) = it^— P<pY ^\ I dxFY N (s,t(s,x)), (10) 
62 its Pkn{s) J 

t(s, x) = ml + m 2 K - 2 + pj y (s) + p 2 KN (s) + 2p KN (s) p y (s) ar , 

P 2 kn(s) = i((s-m%- M 2 N ) 2 — 4 m 2 K M 2 N ) , pJ y ( S ) = i(( a - m 2 - M 2 ) 2 - 4mjM 2 ) . 

The isospin averaged cross section for the if iV — > 0A reaction is equal to 

0- K N-,cf>A{s) = -O-ajv(s) (11) 

where we use that &k-p-->4>a( s ) — a K°n^k{ s ) ■ Evaluating the isospin averaged cross 
section for the KN — > reaction we take into account that the reactions with charged 
E hyperon in the final state have additional isospin factor 2, then 

&KN^M S ) = a XN(s) . (12) 

The cross sections (fTTj) and (fl2"l) are shown on the right panel in figure [2j In comparison 
with the 7r~p — > (pn cross section the KN — > cross sections are about 60 times larger 
at As 1 / 2 ~ 100 MeV. 



3. Strangeness production 



In order to estimate the efficiency of reactions with the (j) production on hyperons 
and kaons, we have to study, first, the dynamics of the strangeness production in 
the course of a nuclear collision. This question was addressed within hadro-chemical 
models in [39-42,55]. In our consideration we make the following assumptions: (a) the 
fireball matter is baryon-dominated; (b) strangeness can be considered perturbatively, 
i.e, the number of produced strange particles, say kaons, is much smaller than the baryon 
number. In the baryon-dominated matter the particles carrying strange quarks, K~ , 
K° , A , S and heavier hyperons have short mean free path. Hence, we can assume 
that they are confined inside the fireball and remain in thermal equilibrium with other 
species till the fireball breaks up. Oppositely, the particles with anti-quark , K + and K° 
mesons, have larger mean free paths and can leave the fireball at some earlier stage of the 
collision. Based on this observation the strangeness separation scenario was formulated 
in refs. [44-46]. It can be suitable for nucleus-nucleus collisions at AGS energies. The 
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analysis in [46] shows that it can be also partially valid at SPS energies. However, at 
what stage of the collision kaons (K + and K°) leave the fireball is not important for 
our further calculations if the concentrations of strange particles remain small. Since in 
such a case the reactions with strangeness annihilation, KY — > ttN and KK — > tttt, can 
always be neglected and the number of strange particles accumulated in the system is 
equal to the number of produced kaons. 

The time dependence of the fireball temperature and baryon density will be 
parameterized in the form of a scaling solution of the hydrodynamic equations [47] 

m = i) Q ' pB{t) = WJW+W^ 2 ' (13) 

where T m and p m are the initial (maximal) temperature and density of the fireball, and 
to is the typical time scale of the fireball expansion. 

The evolution of the kaon density is described by the differential equation 

^-p K (t) P J§ = ll(t), PK (0) = 0. (14) 
dt p B {t) 

The second term on the left-hand side takes into account the reduction of the kaon 
density due to the expansion of the system (dilution). We exploit here the relation 
for the volume change V(t)/V(t) = —psit)/ ps(t) . The kaon production rate on the 
right-hand side is determined by the processes with ttN, it A, NN, inr and iVA in the 
initial states 

K = P7T p N + p. p A + K*N P 2 n + Pi + 4a Pn Pa ■ (15) 



As argued above we do not include on the right hand side in ffl4|) the rates of kaon 
annihilation processes. They are clearly absent if K mesons leave fireball immediately 
after their creation. If K mesons stay in the system the annihilation reactions are still 
very suppressed as being proportional to (pk/pb) 2 <C 1 . 
The transport coefficient 



K 



KX = (Vgb Xv ab) 
"'' ^ l + 5 a 



>ab 

is the cross section a^ b x (s) of the binary reaction a + b — » K + X averaged with the 
particle relative velocity v a b over the momentum distributions of colliding particles. 
Following ref. [39] we write 

/oo 
dxa« x (x 2 ) K x (%) [x 2 - (m a + m b ) 2 } [x 2 - (m a - m b ) 2 ] 

ab 1 ) (l + 5 ab )4mlmlTK 2 (m a /T)K 2 (m b /T) ' 1 ) 

where K\(x) and K 2 {x) are the MacDonald's functions of the first and second orders, and 
y/s th is the reaction threshold. Among the ttN induced reactions we include nN — > K A , 
txN -> KY, and nN -> KKN processes: = + + k^ n . For the isospin 

averaging we use the relations [36,48] 

o(nN — > KK) = -a(7r~p^ AK°) , 
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o-(uN -> UK) = l - (a(n-p -> + -> E°X°) + ^vr+p -> E+X+)) , 

a(nN -> iVM) = 3 (r(7r - p -> NK°K~) . 

The parameterizations for the elementary 7rA^ — > XY cross sections are taken from 
refs. [49,50] and for the 7r~p — > NK°K~ cross section from ref. [36]. In the 7rA channel 
we include only the process with ttA and 7rE in the final state: k^a = ^a + K wA ■ 
The isospin averaged cross section of the nA — > FX reactions we derive at hand of the 
relations (21-31) in ref. [49] 



a(nA ^ EX) = ^ (5 a{^~ A++ -> E°X+) + 6 a(7r°A° -> E" 



7^ 



+ 21 a(7T + A° -> E°X+) + 8 a{ji + /Sr -> E~X + )) 
o"(7rA -> EX) = i (7(71- A++ -> AX) 

3 

with the elementary cross section given in ref. [49]. The transport coefficient of kaon 
production in the NN channel has three contributions k^n = k-nn N + k nn N + k nn NN 
and the parameterisations of the elementary NN — > KNY and iViV — > KKNN cross 
sections are taken from ref. [51]. The cross section for the -kit — > i^K reaction which 
determines we take from ref. [51] too. The cross sections for reactions NA — > KYN 
which contribute to k^* = k na K + k wa E are taken from ref. [52] 

In our model consideration we neglect contributions from other process with kaon 
production which have higher thresholds and, therefore, are less relevant in the range 
of temperatures we consider below T < m n . 

The interaction among pions, nucleons and deltas is strong enough to sustain them 
in thermal and chemical equilibrium during the whole fireball evolution. Hence, their 
densities are given by 

p*(t) = r ^ p 2 ^ {-^ m i+p 2 /T(t)) , 

PnH) = ^l°° dpp 2 exp (pL B (t)/T(t) - ^m 2 N+P yT(t)) , 

PA(t) = ^ jH dpp 2 exp (p B (t)/T(t) - yJml+jfilTitj) . (17) 

We disregard the finite width of As and treat them as stable particles with the mass 
mA = 1232 MeV. Neglecting the contributions of hyperons, heavier resonances and 
anti-particles we can find the baryon chemical potential from the following relation 

^2 



p B (t) = ~T(t) In 



_2_ r°° dpp_ / -y/^^/TQ) +Ae -y/ri£&/nt) 

nit) JO 7T 2 V 



(18) 



In order to study different fireball evolution scenarios we consider three sets of 
parameters for the maximal temperature and density 

I. T m = 110MeV, p m = 4p , a = 0.2, 

II. T m = 130 MeV , p m = 5 p , a = 0.3 , 

III. T m = 150MeV, Pm = 6 Po , a = 0.5 , (19) 
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where po = 0.17 fm~ 3 is the nuclear saturation density. The dependence of the 
temperature and density on time is exemplified on the left panel in figure HI On the 
right panel in figure H] the rates of kaon productions in various processes are shown. The 
most efficient ones are processes with irN and NA in the incoming channel. 

When all ingredients entering the rate (|T5l) are specified the density of produced 
kaon follows after integrating ( fl4f) 

,K{t') tf f\ n{T{t x),p B {t x)) 



PK (t) = PB (t) / dt'^-L = PB (t)t / dx — u ^ (20) 

The produced kaons escape from the fireball leaving behind a medium with some net 
strangeness accumulated in form of K, A and E. (H and Q can be neglected since 
their abundances are of the higher orders in terms of the small parameter Pk/pb-) The 
occupation of particles caring the strange quark is according to their statistical weights 

n(t) 

Piit) =p K {t)- — , i = K,A,E, 

n R {t) = / o °° ^ exp (-y/m* K +p>/T(tj) , 
n A (t) = exp (fi B (t)/T(t) - y/ml+i?/T(tj) , 



n E (t) = 3 / ^PjL exp (p B (t)/T(t) - y/ml+p>/T(tj) . (21) 



7T 

dpp 

71 

We assume, thereby, that the reactions redistributing the strange quark among K and 
hyperons are swift. 

The evolution of strangeness in the course of a collision is illustrated in figure [5] 
where we depict the ratios of strange particles, K , K , A , and E, to the baryon density 

r H (t) = p i (t)/p B (t). (22) 



Catalytic <\> meson production in heavy-ion collisions 



12 



2.0 | — , — i — , — i — , — i — , — i — , — | 5 




0.0 0.2 0.4 0.6 0.8 1.0 0.0 0.2 0.4 0.6 0.8 1.0 0.0 0.2 0.4 0.6 0.8 1.0 



*\ m o fto 

Figure 5. Evolution of the number of kaons (K) produced in the course of the collision 
and the number of strange particles (S, A and K) accumulated in the fireball. Three 
panels (a-c) corresponds to the different sets of the maximal temperature and density 
(I-III) in ((Til). 



Note that in view of eqs. ( |T3i) and ( l2TT) the kaon density ( l20|) and the densities of 
other strange particles are linearly proportional to the parameter i f° r the fixed t/to. 
Presenting the results in figure [5] we take into account this scaling, i.e., the lines are to 
be scaled by factor t /(10 fm). 

We observe that strangeness is dominantly presented in form of £ hyperons and 
that the number of produced strange particles doubles roughly with an increase of the 
temperature by every 20 MeV from case I to case III in (fl9l) . 



4. $ mesons in heavy-ion collisions 

In this section we discuss how the catalytic reactions can influence various aspects of 
production in heavy-ion collisions. 



4-1. $ meson production 

First we consider whether the catalytic reaction can be an efficient source of mesons 
compared to a conventional reactions. The rate equation for the meson production 
has the form similar to (EG 



- p * {t) = R ^ + £ R Ut) + XL® + ■ ■ ■ . 

at p B {t) y,y=A,s 

where on the right hand side we explicitly write down the production rates from the 
conventional ttN — > <pN reactions B% N (t) = k% n p w Pn and the catalytic reactions on 
hyperons R^ Y = Dy,y=a,£ Ryy > -^yy(^) = K ty Pn Py and on anti- kaons Rg N (t) = 
( k ^tv + ^Rn) PR Pn ■ The ellipses stand for further <p production reactions like those 
discussed in [29] and the (f> absorption processes, e.g., (pN — ^ KY and (f>Y — > rrY, etc. 
The rates of various process are shown figure O The production in tcN collisions 
starts of course at the very beginning and gradually falls off as the fireball cools down 
and expands. The rates of catalytic reactions increase, the more strange particles are 
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Figure 6. Production rate of <f> mesons in various catalytic reactions and their sum 
(solid line) in comparison with the irN — > (f>N reaction rate (dashed line). Different 
panels (a-c) show the results for sets I— III, respectively 



produced, reach the maximum at (0.3-0.4) t and drop off later. For set I with the 
maximal temperature T m = 110 MeV the net rate of catalytic reactions yields to that 
of ttN — > <f)N. For set II with T m = 130 MeV the rates become comparable for times 
> 0.6 to- For T m = 150 MeV in set III the catalytic production of exceed the rate of 
the irN — > <pN reaction already at 0.2 to- Note that the numerical values for the rates 
in figure [6] correspond to the fireball expansion time to = 10 fm. If the collision lasts 
longer then the curves for catalytic reactions have to be scaled up by a factor t / (10 fm), 
since the number of the strange particles is proportional to the expansion time. This 
will make the catalytic reaction efficient even for the smaller temperature in set I. 

The above estimates show that the catalytic processes can be potentially important 
for the production in heavy-ion collisions. It is interesting now to investigate how they 
can affect other <fr production patterns, and whether their dominance can lead to any 
contradiction with experimental data. 

4-2. Centrality dependence 

We discuss now the centrality dependence of the production. Changing centrality 
one changes the volume of the system. We will use the mean number of projectile 
participants, N pp , as the measure for initial volume of the fireball created in the collision, 
V oc N pp . The N pp can be directly related to the energy deposited in the zero-degree 
calorimeter, cf. [2]. If there is only one changing parameter with the unit of length 
as in the case of a symmetrical collision at a fixed collision energy, I ~ V 1 ^ 3 oc N^ 3 , 
the scaling properties of hydrodynamics imply that the collision time is of the order 
to ~ l/c oc N^ 3 [55]. The multiplicities of particle species which reach the full chemical 
equilibrium, e.g. pions and As, are proportional to the volume of the fireball and scale, 
therefore, as N n oc V oc N pp and ~ iV pp . The multiplicity of kaons (K + ), which 
have a long mean free path and leave the system right after being produced, scales as 
N K + oc V t oc iVp 1 / 3 . According to (|2T]) other strange particles should follow the K + 
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Figure 7. The <f>/ir (left panel) and 4>/K + (right pane) as a function of the mean 
number of projectile participants. The data points from ref. [2] correspond to Au+Au 
collisions at 11.7 AGeV. The curves show the relations (|2U)) . The value of the 
parameters are given on the plots (A = 197). 



scaling, N^^k ~ ^pp 3 • The number of produced mesons can be estimated as 



N n N p „ N^Ny . N R N PP 



AT" " PP ± I 71 2 j . _ ~ ' /\ ~ ' pp , 

iV^ ~ K^jv y r + K 7ry y r + K KN y Z 



~ ^conv -^pp "I - ^cat -^pp 

The term ~ ^pp 3 is due to the conventional production reactions like tcN — > cf)N whereas 
the term ~ one corresponds to the catalytic reactions which contribute with an 
additional factor t . For the ratios experimentally observed in ref. [2] we find 

where a , a' and b , b' parameterize the relative strength of conventional and catalytic 
processes, and A is the number of nucleons in the colliding nuclei. In figure [7] we 
compare these parameterizations with the only available data of the yield centrality 
dependence for Au+Au collisions at 11.7 AGeV [2]. First, we adjust parameters a, a' 
and b, b' separately and obtain dashed and solid curves, respectively. Because of large 
error bars both fits with a, a' = or b, b' = equally well go through data points. So 
the data on centrality dependence cannot rule out even the dominance of the catalytic 
reactions. The optimal fits are reached when both parameters are activated, dashed- 
dotted lines in figure [71 Comparing the values of parameters a, a' and b, b' from these 
fits we can conclude that the catalytic reaction contribution can be about 30%-40% for 
iV pp = A. 



4-3. (j) rapidity distribution 

The systematics of rapidity distributions in Pb+Pb collisions at the SPS is reported 
in ref. [7]. The distributions can be fitted with a sum of two Gaussian functions placed 
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Figure 8. Root mean square of the rapidity distributions of particles produced in 
Pb+Pb collisions at SPS energies versus the beam rapidity. Left panel: experimental 
data for tt~ and K from ref. [53,54], <ft mesons from [7] and A from [56]. Right 
panel: The distribution width for <fi mesons produced in 7rA — > 4>Y, K + A — > (f>N and 
K + K ~ — > <f) reactions calculated with eq. (flZu) , except the dash-dot-dot line calculated 
with eq. (|25|) in comparison with the experimental data from [7] . 



symmetrically around mid-rapidity 

1 diV 1 r (»-°) 2 (y+g) 2 1 , , 

(JV> dy V8tt(t 2 L J 
here y is the rapidity in cm system. The width of the distribution is characterized 
by the root mean square RMS 2 = a 2 + a? . The left panel of figure [H] depicts the 
RMSs versus the beam rapidity for mesons, constructed from the parameters a and 
a given in table V in [7], together with those for ir~ , K + and K~ [53,54]. Ref. [7] 
pointed out that the width of the 4> meson distribution does not fit into the systematics 
observed for the other mesons but increases much faster with the collision energy. It was 
emphasized in [7] that the steep rise of the <fi distribution width cannot be explained by 
the hadronic process K + K~ — > <p. The rapidity distribution of from the latter process 
can be estimated as the product of K + and K~ rapidity distributions. The width of 
the product was calculated in [7] under the assumption of the single-Gaussian sources 
for kaons with RMSx+ and RMS^-- as 

RMS^-= RMSg+ RMSff- 
yRMS K + + RMS^- 

We show this result on the right panel of figure [8] by the dash-dot-dot line. The 
resulting width is much smaller than that observed in the experiment. Question is 
if any other <fr production mechanism can produce the steep rise of the RMS with the 
collision energy. We note that the rapidity distribution of hyperons increases much 
faster than those for mesons as the hyperons are dragged with nucleons to forward and 
backward rapidities (an effect of partial transparency). The RMSs of A reported in 
ref. [56] are shown on the left panel in figure M by solid circles. The distribution width 
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for 158 AGeV was not determined experimentally, therefore we calculated it from the 
rapidity distribution obtained theoretically in [57] within a hydrodanamic model which 
successfully reproduces the particle production. 

We assume that the rapidity distributions of particles do not change after some 
initial stage when nuclei are passing throuhgh each other. This implies the absence or 
weakness of acceleration and diffusion processes. The collision kinematics is restricted 
mainly to the exchange of transverse momenta. Then the rapidity distribution of 0s 
produced in the reaction 1 + 2— > 0+ AT is roughly proportional to the product of rapidity 
distributions of colliding particle species 1 and 2. If the distributions are given by the 
two Gaussian distributions (12411 the RMS of the resulting <p distribution is given by 



RMSi 2 



dyy 2 / (y, a 1; oi) f(y,a 2 ,a 2 ) 



dyf{y,a h ai) f(y,a 2 ,a 2 ) 



a\a\ + a\a\ + a\o\ + 2 a, a 2 o\o\ ^ (2g) 



This expression differs from (125]) applied in [7]. Using the parameters of K + and K~ 
from [53, 54] and eq. (I26p we obtain RMSs for the K + K~ — > <fi reaction shown in 
figure [8] (right panel) by dotted line. They are somewhat larger than those obtained 
in [7] but still much smaller than empirical ones. In contrast, the width of <p rapidity 
distributions from the reactions involving A particles, 7rA — > <f)Y and K + A — > <pN (solid 
and dashed lines in figure [8]) rises much faster and is comparable with the experimental 
results for beam energies between 20 AGeV and 80 AGeV. In collisions at 158 AGeV 
the <p distribution is much broader than our estimates. Perhaps, some new non-hadronic 
mechanism of the production becomes operative at this energy. 



5. Conclusions 



In this paper we study a new mechanism of meson production in nucleus-nucleus 
collisions — the catalytic reactions on strange particles, e.g., ttY — > cj)Y and KN (pN. 
These reactions are OZI-allowed and their cross section can be by an order of magnitude 
larger than the cross section of conventional OZI-suppressed <fi production reactions 
ttN — > <pN and A^A — > <j)N, etc, considered so far. These reactions require only 
one strange particle in the entrance channel and therefore are less suppressed than 
the strangeness coalescence reactions, KK — > <p and KY — > <pN in collisions where 
strangeness is statistically suppressed. Using a hadronic Lagrangian ([3]) we estimate the 
nY — > cf)Y and A^/V — > <j)N cross sections to be roughly of the order of 1 mb. In order to 
estimate the efficiency of the new reactions we calculate the evolution of the strangeness 
content of the fireball within a hadrochemical model. The main assumption we follow 
is that kaons (K + and A' ) can leave the fireball freely after being created whereas the 
strange particles (hyperons, K~ and K°) remain in the fireball in thermal and chemical 
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equilibrium. Thus the net strangeness accumulated in the fireball is growing in the 
course of the collision. We compare the rates of production in the catalytic reactions 
and in the ttN — > <f>N reaction. The former can be competitive in collisions with the 
maximal temperature above 110 MeV and the collision time > 10 fm. The efficiency 
of catalytic reaction increases if some strangeness is presented in the fireball right in 
the beginning after the first most violent nucleon-nucleon collisions and if the fireball 
lifetime is longer. 

We discuss how the catalytic reaction could affect the centrality dependence of 
the yield. Since the catalytic rates depend on the concentration of strange particles 
accumulated during some time, the resulting number of 0s produced in such a reaction 
would grow oc ig, where tg is the collision time, in contrast to the number of 0s from 
conventional reactions irN — > <pN growing oc to • in the scaling regime of hydrodynamics 
the collision time can be related to the typical length scale of the system determined 
by its volume to oc I oc V 1 ^ 3 oc A^ 3 , where N pp is the centrality criterion - the mean 
number of projectile participants. Taking two types of processes into account we can 
parameterize the experimental centrality dependence of the N^/N^ and N t j ) /N K + ratios 
measured in Au+Au collision at 11.7 AGeV as follows 



Here the first terms are due to the conventional reactions involving non-strange particles 
and the second terms are due to the catalytic production on strange particles. Relative 
strength of two terms shows that the contribution of the catalytic reactions can be up 
to 30%-40% at the AGS energies for A pp = A. 

Analyzing the rapidity distributions at SPS energies we find that the strong rise 
of the distribution width with the collision energies between 20 AGeV and 80 AGeV 
can be explained by the production in nA and K + A collisions. 

The present analysis of the catalytic mechanisms of the production is only 
exploratory and should merely serve as an invitation for further investigations. 
Particularly, a more reliable calculation of the cross sections for the ttY — > <pN and 
KN — > (f>N reactions and fireball evolution is mandatory. 
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